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We analyze the relationship of generalized conditional symmetries of evolution equations to the 
formal compatibility and passivity of systems of differential equations as well as to systems of vector 
fields in involution. Earlier results on the connection between generalized conditional invariance 
and generalized reduction of evolution equations are revisited. This leads to a no-go theorem on 
determining equations for operators of generalized conditional symmetry. It is also shown that up 
to certain equivalences there exists a one-to-one correspondence between generalized conditional 
symmetries of an evolution equation and parametric families of its solutions. 



^ ; 1 Introduction 

Ph. 

Generalized conditional symmetries provide an effective method for finding exact solutions of 
evolution equations. Similarly to other such methods [H], it can be viewed as an instance of 
the general method of differential constraints |60t [65] (or "side conditions" [Sj). Within the 
framework of empiric compatibility theory, generalized conditional symmetries as differential 
| constraints compatible with an initial equation were investigated by Olver [13] in order to justify 

£> ' the method on "nonlinear separation" of variables by Galaktionov [20J . Another interpretation 

of generalized conditional symmetries of an evolution equation is to consider them as invariant 
\ manifolds of this equation, i.e., manifolds in appropriate jet spaces that are invariant under 

the flow generated by the equation. This is the terminology in which generalized conditional 
symmetries of systems of evolution equations were first studied by Kaptsov [IJ [32] although the 
importance of invariant manifolds of evolution equations was understood much earlier [36j . 

From the symmetry point of view, the notion of generalized conditional symmetry arises by 
merging the notions of generalized and conditional symmetries, cf. Section [2j The idea of signifi- 
cantly extending Lie symmetries of differential equations by including derivatives of the relevant 
dependent variables in the coefficients of the associated infinitesimal generators first appeared 
in the fundamental paper of Noether [40J in connection with her study of conservation laws. 
Symmetries of this kind are called, e.g., generalized [12], Lie-Backlund [El EI] or higher-order [6] 
symmetries in the literature. See additionally the excellent sketch on the history of generalized 
symmetries and relevant terminology in [421 p. 374-377]. The concept of conditional symmetries 
arose much later. Its origin can be traced back to the thesis of Bluman [5] and the paper by 
Bluman and Cole [7], where it was presented in terms of "nonclassical groups" or the "non- 
classical" method of finding similarity solutions, respectively, cf. the detailed discussion in [61 
Section 5.2.2]. A version of the corresponding invariance criterion explicitly taking into account 
the differential consequences involved in the process was first proposed by Fushchych and Tsyfra 
in [18]. Combining results of [HI [18] and other previous papers, in [13] Fushchych introduced 
the general concept of conditional invariance. Around this time the terms "conditional invari- 
ance" and "Q-conditional invariance" began to be used regularly in connection with the method 
of Bluman and Cole and soon evolved into the terms Q-conditional [16] or, simply, conditional 
|19j and nonclassical [37] symmetry. The notions of generalized and conditional symmetries 
were merged, within the framework of symmetry analysis of differential equations, by Fokas and 
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Liu [12] in the special case when evolution equations and symmetries do not explicitly involve 
the time variable and by Zhdanov [67] in the general case. 

The variety of possible interpretations and related notions and a number of different names 
for the parent notions of conditional and generalized symmetries leads to the diversity of names 
used for generalized conditional symmetry in the literature. We have already mentioned the 
terms "invariant manifold" [HEJ|32] (resp. "invariant set" |2H I24j) and "compatible differential 
constraint" [43]. Additionally, combining names of the parent notions of symmetries leads, in 
particular, to the terms "conditional Lie-Backlund symmetry" [291 130} 167] and "higher (or higher 
order) conditional symmetry" [HEH]. Sometimes special names are used for particular cases of 
generalized conditional symmetries. For example, linear compatible differential constraints for 
diffusion-reaction equations were called "additional generating conditions" in [9] . For uniformity, 
we will use the term "generalized conditional symmetry" |12U53j, l54j throughout the paper. This 
will additionally emphasize the relation of this notion to symmetry analysis although the nature 
even of usual conditional symmetries is in fact closer to compatibility theory, cf . [35] . 

The main purpose of this paper is to investigate basic problems concerning generalized con- 
ditional symmetry of (1 + l)-dimensional evolution equations of the general form 

ut = H(t,x,U( r>x )), (1) 

where r ^ 1, Ut = du/dt, uq := u, Uk = d k u/dx k , it( r)X ) = (uo,u\, . . . ,u r ) and H Ur ^ 0. Among 
these problems are the comparative analysis of different versions of the conditional invariance 
criterion, the study of the possibility of solving the corresponding determining equations as 
well as relating generalized conditional symmetries to the concept of reduction, multiparametric 
families of solutions and different notions of compatibility for overdetermined systems of partial 
differential equations. Most results of the paper can be extended to systems of (l+l)-dimensional 
evolution equations if certain restrictions for generalized conditional symmetries are imposed, 
cf. p]. We restrict our consideration to single evolution equations for the sake of clarity of 
presentation. 

Throughout the paper we denote by E a fixed equation of the form ([I]) . The indices a and b run 
from 1 to p, and we use the summation convention for repeated indices. Bar over a letter denotes 
a tuple of p consecutive values. Subscripts of functions denote differentiation with respect to 
the corresponding variables, dt = d/dt, d x = d/dx, d u = d/du and u t k = d k+1 u/dtdx k . We also 
will use another notation for derivatives: u a = u a0tCei = d^u/dt a °dx ai , where a = (ao, a\) is a 
multiindex, ao,ai G NU {0} and \a\ = «o + so that Uk = «o,fc an d = ui £. Any function 
is considered as its zero-order derivative. D t = dt + u aQ+ i yai d Ua and D x = d x + d Ua are 

the operators of total differentiation with respect to the variables t and x, respectively. All our 
considerations are carried out in the local setting. 

In the next section we discuss prerequisites for introducing the notion of generalized con- 
ditional symmetries in symmetry analysis and present different versions of the corresponding 
invariance criterion for single evolution equations. Relations of generalized conditional symme- 
tries to formal compatibility and passivity of certain overdetermined systems of partial differ- 
ential equations as well as to involutivity of certain systems of vector fields are established in 
Sections [31 H] and [5l respectively. For this purpose we employ a weight of derivatives instead of 
the usual order (Section [3]) and a ranking of derivatives (Section [4]), which are associated with 
evolution equations of a fixed order. Reductions of evolution equations with special ansatzes 
are studied in Section [6] The Zhdanov theorem [67\ [68] (see also [1]) on the connection of 
generalized conditional symmetries of an evolution equation with ansatzes of a special form re- 
ducing this equation is also revisited. This leads to new results on the correspondence between 
generalized conditional symmetries, ansatzes and parametric families of solutions of evolution 
equations. In Section [7J we prove a no-go theorem on determining equations for generalized 
conditional symmetries of evolution equations. Roughly speaking, it is shown that solving the 
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determining equations is equivalent to solving the original equations. An interpretation of usual 
conditional symmetries of evolution equations as special generalized conditional symmetries is 
given in Section [8] and is then illustrated by a new nontrivial example. Properties of generalized 
conditional symmetries of evolution equations are summed up in the conclusion. 

2 Different forms of the criterion of conditional invariance 

The criterion of generalized conditional invariance of evolution equations arises as a natural ex- 
tension of both the criterion of generalized invariance and the criterion of conditional invariance. 
This is why we at first analyze the latter criteria in the case of evolution equations. 

By the conventional definition, an equation £ of the form ([1]) is conditionally invariant with 
respect to the vector field Q = rdt + id x + nd u , where the coefficients r, £ and rj are functions of 
t, x and u, if the relation Q^E\ £ n g = holds. Here E := ut — H and the symbol Q( r ) stands 
for the standard rth prolongation of the operator Q [42} H6] : 

Q(r)=Q+ E (D?°D^Q{u]+TU ao+ljai +t; Uao>ai+1 )d Uoi , (2) 

0<\a\^r 

where Q[u] = r\ — rut — £,u x is the characteristic of the vector field Q, and Q r denotes the 
manifold defined by the set of all the differential consequences of the characteristic equation Q: 
Q[u] = in the rth-order jet space J r , i.e., 

Q r = {(t,x,u {r) ) G J r | D£°D% l Q[u] =0, Q + ai < r}. 

The manifold defined by the equation £ in J r is denoted by £ r . In comparison with classical 
Lie symmetries, the weakening of the invariance condition consists in equating Q^E to zero on 
the submanifold £ r n Q r but not on the entire manifold £ r . As £ is an evolution equation, only 
differential consequences of Q with respect to x are in fact essential when substituting into the 
expression Q/ r \E (cf. the proof of Proposition Hence the conditional invariance criterion can 

be rewritten in the form Qr r \E\ c _ =0, where 

it)-ny( ria .) 

Q(r,x) = {(t,x,u (r) ) E J r | D k x Q[u] = 0, k = 0,...,r- 1}, 

and the bound r for orders of the occurring differential consequences of the equations £ and Q 
is not essential. 

Two vector fields Q and Q are called equivalent if they differ by a multiplier which is a 
nonvanishing function of x and u: Q = XQ, where A = X(x,u), A ^ 0. The property of 
conditional invariance matches nicely with this equivalence relation. Namely, if the equation £ 
is conditionally invariant with respect to the vector field Q then it is conditionally invariant with 
respect to any operator which is equivalent to Q. Therefore the equivalence relation of vector 
fields has a well-defined restriction to the set of conditional symmetries of the equation £ . 

In the case of generalized symmetries, the extension of the notion of Lie symmetries is to 
permit the dependence of coefficients of vector fields on derivatives of u |42j . A generalized 
vector field Q is a symmetry of £ if and only if the associated evolutionary vector field Q[u]d u 
is. Hence it is sufficient to consider only evolutionary vector fields as generalized infinitesimal 
symmetries. Additionally, if an evolutionary vector field Q = rjd u is a symmetry of £ and the 
difference fj — n vanishes on solutions of £ then the vector field Q = rjd u also is a symmetry of £. 
Such generalized symmetries are called equivalent. In view of the evolution form of £ this means 
that we need to consider only generalized symmetries whose characteristics do not depend on 
derivatives containing differentiation with respect to t. 

Merging the above extensions of classical Lie symmetries leads to the notion of generalized 
conditional symmetries. Consider a generalized vector field Q = rjd u with r\ = Q[u] being a 
differential function, i.e., a smooth function of i, x and a finite number of derivatives of u. 
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Definition 1. An evolution equation £ is called conditionally invariant with respect to the 
generalized vector field Q = r]d u if the condition 

Q{r) E \ M = 

holds, where the rth prolongation Qr r \ of Q is defined by ([2]) and M. denotes the set of all 
differential consequences of the equation £ and differential consequences of the equation n = 
with respect to x. In this case, Q is called an operator of generalized conditional symmetry of 
the equation £ and the above condition is the criterion of conditional invariance. 

As Qi r )E = Dtt] — Ylk=o ^u k D^.r] and the last sum identically vanishes in view of differential 
consequences of the equation n = with respect to x, we obtain at once another form of the 
criterion of conditional invariance |67j : 



After calculating the orders of the occurring differential consequences, Definition Q] can be 
equivalently reformulated with a precise determination of the underlying jet space. To this 
end, it suffices to consider the criterion Dti]\ M = within the jet space J rn of order 
m = max{r(ao + 1) + «i | r\ Ua ^ 0} which coincides with the weight of Dffj (cf. Section[3]). Then 
the criterion takes the form Dtn\ Mm = 0, where Ai m is the manifold determined by Ai in J m . 
All other similar conditions can be formalized in the same way. 

There are two well-defined equivalence relations on the set of generalized conditional sym- 
metries of the equation £, which extend the above equivalence relations of conditional and 
generalized symmetries, respectively. 

Suppose that fj = Xn, where A is a nonvanishing differential function, i.e., Q = nd u and 
Q = fjd u are equivalent generalized vector fields. Then the vector field Q = rjd u is a generalized 
conditional symmetry of the equation £ if and only if the vector field Q = fjd u is. Indeed, 
Difj = XDtr] + r/DfX vanishes assuming Ai if and only if Dtrj does. Moreover, Dtfj vanishes 
assuming Ai if and only if it vanishes assuming Ai, where Ai denotes the set of all differential 
consequences of the equation £ and differential consequences of the equation fj = with respect 
to x. This allows one to restrict the equivalence relation of generalized vector fields to the set of 
generalized conditional symmetries of the equation £ in a well-defined way, analogously to the 
case for usual conditional symmetries. Hence we will say that generalized conditional symmetries 
Q = V^u and Q = fjd u of £ are equivalent as vector fields if there exists a nonvanishing differential 
function A such that fj = Ar/. 

If differential functions r/ and fj coincide on the manifold defined by differential consequences 
of £ in a jet space of suitable order, then in view of the Hadamard lemma we have a representation 
fj = 7] + x a Df°D^ 1 E, where the summation is over a finite set of a's and the x a are differential 
functions. Hence the condition D t r]\ M = is equivalent to the condition D t fj\ M = and, 
therefore, the condition Dtfj\j^ = 0. In other words, the vector field Q = nd u is a generalized 
conditional symmetry of £ if and only if the vector field Q = fjd u is. For this reason we will call 
the generalized conditional symmetries Q = nd u and Q = fjd u equivalent on solutions of £. 

In contrast to the equivalence of generalized conditional symmetries as vector fields, the equiv- 
alence on solutions does not agree with the reduction procedure. Some vector fields from a set 
of generalized conditional symmetries equivalent on solutions of £ cannot be used for reduction 
of £, while some of them are appropriate for reduction but the corresponding reduction proce- 
dures differ in the number of invariant independent and dependent variables in the associated 
ansatzes and, therefore, the structure of the reduced systems, cf. Section [8l 

We can merge the above two equivalence relations of generalized conditional symmetries into 
a single notion. Namely, generalized conditional symmetries Q = r]d u and Q = fjd u of £ are 
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called equivalent if there exists a nonvanishing differential function A such that fj — Xrj is equal 
to zero on solutions of £. 

Taking into account the equivalence on solutions of the evolution equation £ , we can restrict 
our considerations to generalized conditional symmetries of the reduced form Q = fjd u , where 
the characteristic fj is a reduced differential function, i.e., it depends on t, x and derivatives of u 
with respect to only x. Generalized conditional symmetries in reduced form are equivalent if 
and only if their characteristics differ in a nonvanishing multiplier being a reduced differential 
function. Up to this equivalence, we can replace Q by the corresponding canonical form 

Q = {u p - f)(t, x, u^ liX) ))d u , (3) 

where p is the order of fj and the condition of maximal rank of fj with respect to u p is additionally 
assumed to be satisfied. The function fj = v{t,x,u^ p _ lx ^) is obtained by solving the equation 
fj = with respect to u p . 

An evolution equation £ is conditionally invariant with respect to a generalized evolution 
vector field Q = rj(t, x, ur PjX ))d u in reduced form if 

Q{r){ut-H)\ =0, or D t r]\ =0, (4) 

where Qt r ) is the rth prolongation of Q defined by (|2|), £ r + P (resp. Q( r +p,x)) is the manifold 
determined in the (r + /j)th-order jet space by differential consequences of the equation £ (resp. 
the equation r\ = only with respect to x). HQ is in canonical form, i.e. ij = u p — ij(t, x, Uf p _i iX \), 
the criterion of conditional invariance of £ with respect to Q reads 

D p x H = D t fj on {u p+k = D*fi, k = 0,...,r, u tl = D l x H, I = 0, . . . , p - I}. (5) 

After making all necessary substitutions in ([5]), we obtain the single determining equation 

DtV = D p x H (6) 

in fj, where H = H(t, x, no, . . . , u r ) if p > r, H = H(t, x, no, • • • , u p -i,f], D x f], . . . , D r x p fj) if 
p ^ r, and 

D t = d t + (D b x - 1 H)d Ub _ 1 , b x = d x + Y, u b d Ub ^ + fjdu^ 

6=1 

are the operators of total differentiation restricted to the manifold £ r + p r)Q( r+p >x y Equation (|6|) is 
a (1 + /^-dimensional evolution equation in an unknown function fj of the independent variables 
t, x, uq, . . . , u p -i, and we have no possibilities for splitting with respect to unconstrained 
variables. 

There also exist other forms and interpretations of the criterion of generalized conditional 
invariance of evolution equations in the literature. Suppose that the generalized evolution vector 
field Q is in reduced form. On the manifold £ r + P we have Qr r )(ut—H) = vt+V*H— H^r], where /* 
denotes the Frechet derivative of a differential function / depending solely on t, x and derivatives 
of u with respect to x, 



i=0 



Since the differential function rjt + v*H — H*r] does not involve derivatives with respect to t and 
mixed derivatives, we can rewrite in the form 

(Vt + V*H - H*-n)\ =0, or (r? t + v*H)\ Q =0. 
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If rj t = rj x = 0, H t = H x = and n is of maximal rank with respect to u p , in view of the 
Hadamard lemma the last condition is equivalent to the condition rj*H — H*rj = F[u, rj\ presented 
in Definition 1.1 of [12] . Here F[u,rj] is a smooth function of derivatives of u with respect to x 
and total derivatives of rj with respect to x such that F[u, 0] = 0. 

Introducing the notation Dt for the reduced operator of total differentiation with respect to t 
on the solution set of the equation £, 



we represent T]t+r)*H as Dtr\ and obtain as another form of the criterion of generalized conditional 
invariance of evolution equations 



which can be interpreted as the condition of invariance of the equation rj = with respect to the 
formal transformation group [27] generated by the generalized vector field Dt. Since the vector 
field D t is associated with the equation £, the solution set of the equation rj = is called an 
invariant set, or, interpreted as a manifold in an appropriate jet space, an invariant manifold 
of the equation £ [TJ Section 3.1]. This interpretation is especially clear in the case rjt = and 
Ht = 0. Then we can rewrite the criterion in the form 



consider t as the group parameter of the formal transformation group corresponding to the 
generalized vector field Hd u in evolution form and interpret the equation £ as the equation for 
finding this group. 

Remark. Both symmetries and cosymmetries of an evolution equation are generalized condi- 
tional symmetries thereof but they obviously do not exhaust the entire set of its generalized 
conditional symmetries. For example, countable sets of independent symmetries and conserva- 
tion laws had been known for the Sawada-Kotera equation ut = u§ — 30uu^ — 30uiU2 + 180« 2 «i. 
Recently a series of generalized conditional symmetries of this equation, which are neither sym- 
metries nor cosymmetries, was explicitly constructed in [3]. 

3 Formal compatibility and conditional symmetry 

The relations between usual conditional (nonclassical) symmetries, reduction and compatibility 
of the combined system consisting of the initial equation and the corresponding invariant surface 
equation were discovered in [50] and were also studied and extended to the generalized framework 
in [33]. In particular, it was shown that the conditional invariance criterion is the compatibility 
condition of the combined system. This also was remarked, e.g., in [12) . At the same time, 
the rigorous formalization of this relation is nontrivial and was not considered so far even for 
evolution equations. 

In this section we use the definition of formal compatibility as presented, e.g., in |47 [ 1561 [57] . 
We temporarily employ notations compatible with these references, hence slightly different from 
the rest of the paper. 

Let Ck be a system of I differential equations L l [u] = 0, . . . , L l [u] = in n independent 
variables x = {x\, . . . , x n ) and m dependent variables u = (u , . . . , u m ), which involves deriva- 
tives of u up to order k. The system is interpreted as a system of algebraic equations in 
the jet space J k and defines a manifold in J k , which is also denoted by The sth-order 
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D t = d t + J2(D k x H)d Uk , 



k=0 
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prolongation C^+s of the system Ck, s G N, is the system in J fc+S consisting of the equations 
D" 1 . . . D% n LP[u] = 0, j = 1, . . . ,1, \a\ ^ s. Here D{ is the total derivative operator with respect 
to the variable X{. The projection of the corresponding manifold on J k+S ~ q , where q G N and 
q ^ s, is denoted by C^ s _ . The system £fc is called formally compatible (or formally integrable) 
if 4+ s = £fc+« for any s G N U {0} 071 [561 [57]. 

The first obstacle in harmonizing the above definition of formal compatibility and the defini- 
tion of generalized conditional symmetry of evolution equations is that the equations £ and r\ = 
have, as a rule, different orders. Therefore, trivial differential consequences of these equations 
should be attached to the joint system of £ and r] = before testing its compatibility. 

The other obstacle is that the order of each of these equation may be lowered on the manifold 
of the other equation. To avoid this, we take the following steps. 

Firstly, we replace the equation rj = by the equation f] = which is equivalent to the 
equation 77 = under the condition that £ is satisfied, does not contain derivatives involving 
differentiation with respect to t and is of minimal order among equations possessing these prop- 
erties. In other words, we convert the generalized vector field Q = rjd u into its reduced form 
Q = V^u, where fj is of minimal order. 

Secondly, instead of the usual order of derivatives and differential functions with the inde- 
pendent variables t and x we use the weight w defined by the rule: 

w(t) = w(x) = 0, w(u a ) = [a] := rao + a±. 

The technique of working with a weight does not differ essentially from the order technique and 
so a number of analogous notions can be introduced. Thus, in the weighted jet space J^(t,x\u) 
we include the variables whose weight is not greater than k. The weight w(L) of any differential 
function L = L[u] equals the maximal weight of variables explicitly appearing in L. The weight 
of the equation L[u] = equals w(L). In particular, w(ut) = w(H) = r. This implies that the 
weight of the equation £ cannot be lowered by using differential consequences of the equation 
fj = 0. The introduction of the weight also justifies the exclusion of the derivative Ut and mixed 
derivatives from 77 since in contrast to the usual order the weight cannot be raised under this 
exclusion. Note that the weight is also preserved by admissible transformations of evolution 
equations. As for any point or contact transformation between two evolution equations the 
expression of the transformed t depends only on t [331 [38] , and the weight of every differential 
function L[u] is invariant with respect to such transformations. 

Given a system Ck of I differential equations L 3 [u\ =0, j = I,..., I, in the independent 
variables (t, x) and the dependent variable u, which involves derivatives of u up to weight k, the 
sth weight prolongation Ck+ S of the system Ck, s G N, is the system in J k+S (t,x\u) consisting 
of the equations Df° D^ 1 LP [u] = 0, [a] ^ s. The system Ck+s is constructed from the system 
Ck+ s -i by attaching to Ck+ s -i the equations D^°D^. l U[u] =0, [a] = s. The set of these 
attached equations can be viewed to consist of the equations obtained via acting by D x on 
Df° D" 1 LP : [u] = 0, [a] = s — 1, and, if r divides s, the equation obtained from D s J r 1 LP[u] = 
via acting by D t . 

Let s = max(r, p), i.e., s is the weight of the joint system S of the differential equations £ 
and i) = 0, where p = w(fj) = ordr). Denote by P q and V q , where s, the system 

D*fj = 0, k = 0, . . . ,q- p, Df°D^{u t - H) = 0, [a] < q- r 

of algebraic equations in the jet space J$(t, x\u) and the corresponding manifold, respectively. 
In particular, the system P s is obtained via completing the reduced systems of £ and fj = by 
trivial differential consequences which have, as equations, weights not greater than s. 

Proposition 1. The system P s is formally compatible if and only if the evolution equation £ is 
conditionally invariant with respect to the operator Q = rjd u . 
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Proof. By R q , where q ^ s, we denote the (q — s)th weight prolongation of the system P s . Thus, 
the system R s coincides with P s . Additionally to the equations of P q , the system R q includes 
the equations DfD^fj = 0, where [a] ^ q — p and oiq ^ 0. 

Suppose that the system R s is formally compatible. Consider the differential function 

F = D t fj- H Ur D r x f, - f, Up DP(u t - H). 

The equation F = is a consequence of R r + P , and w(F) ^ r + p — 1. As -fi^+p-i = Rr+p-l 
by assumption, the equation F = also is a consequence of the system i? r+p _i which coincides 
with the system P r+p -\. We conclude that F\-p r+p _ 1 = and, therefore, D t f)\-p r+ = 0. The last 
equality is nothing but a form of the conditional invariance criterion. 

Conversely, let the evolution equation £ be conditionally invariant with respect to the operator 
Q = r/du- Then we prove by induction with respect to the value q that R q = P q . The equality 
is obvious for q = s. Supposing that the equality is true for a fixed q, let us prove it for 
q + 1. As R q = P q , the prolonged system R q +\ includes P q +i as a subsystem and additionally 
contains the equations DtD x fj = 0, I = 0, . . . , q + 1 — p — r, which are identities on V q +\ since 
DtD l x fj\-p r+ +l = D l x Dtfj\"p r+ +l = 0. (To prove this last equality, use the fact that D t f]\-p r+p = 0, 
apply the Hadamard lemma, and act by D l x on the resulting representation.) Hence R q +i = Pq+i, 
completing the induction. Among the left hand sides of equations from P q +\ only the differential 
functions D x +1 ~ p fi and D^D" 1 (ut — H), [a] = q + l — r depend on variables of weight q + 1, and 
they are functionally independent with respect to these variables. Hence R<p = P q = R q . □ 

4 Passivity and conditional symmetry 

For the convenience of the reader, at first we briefly present basic notions of Riquier's compat- 
ibility theory. See, e.g., [39] and references therein for a more extended presentation of these 
notions and related results. We again use the notation from the beginning of the previous sec- 
tion. In what follows the indices a and b run from 1 to to, the indices i and j run from 1 to n, 
a and (3 run through the multiindex set {(a\, . . . , a n ) \ oti G N U {0}}. 

Usually the set of derivatives {u%\ is assumed partially ordered. A derivative u a a is said to 
be lower (resp. strictly lower) than a derivative u^, and we write u a a ^ u b ^ (resp. u% < w0, 
if a = b and «j ^ ft (resp. a = b, a% ^ ft and a ^ (3). In contrast to this, the initial 
point of Riquier's theory is a suitable total ordering of derivatives, which is compatible with 
differentiations. Namely, a ranking is a total (or linear) ordering =^ of derivatives such that 
!i„ -< DiU a a and if u a a -< u b p then DiU a a -< DiiA. (As usual, u a a -< means that u a a =<! and 
u a 7^ u /3-) I n view of these properties of a ranking, the condition u a a ^ implies u a a =<! u». 

Suppose that a ranking of derivatives is fixed. By the leading derivative of a differential 
function F[u] we mean the maximal element in the finite set of derivatives {u^ \ F u a ^ 0} if 
this set is not empty. Consider a system C of finitely many differential equations resolved with 
respect to their leading derivatives: 

v% B = F'[u], 8 = 1,..., J. 

The set of leading derivatives of C consists of the leading derivatives of the above equations, i.e., 
it equals {u^ \ 3uJ s : = n^}. The infinite prolongation of the system C is formed by all 
possible differential consequences 

Each of the differential consequences is automatically resolved with respect to its leading deriva- 
tive, which is called a principal derivative of the initial system C. In other words, the set of 
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principal derivatives of C consists of the derivatives of the leading derivatives of C The other 
derivatives are called parametric derivatives of C. 

Differential consequences of C involving only parametric derivatives are said to be integrability 
(or compatibility) conditions. A system C is active if it has unsatisfied integrability conditions, 
otherwise it is called a passive system. 

A system C of equations resolved with respect to its leading derivatives is called 

• triangular if every leading derivative of C is the leading derivative of only one equation 

• autoreduced if no principal derivative occurs on the right hand side of any equation of C 

• orthonomic if it is triangular and autoreduced. 

It is obvious that all of the above properties depend on the choice of ranking. 

Let us return to evolution equations of the form ([1]) . The basic idea for introducing a ranking 
is to assume that u r -< ut -< u r +i- The extension of the last condition to all derivatives of u 
leads to the following ranking: 

u a 4up <^=> [a] < [/3] V ([a] = [0] A a < A))- 

We recall that [a] = rao + a±. This ranking agrees well with the derivative weight introduced in 
the previous section. We rank derivatives by their weight and then use the lexicographic order 
for derivatives with the same weight. 

After this ranking is fixed, the exclusion of derivatives involving differentiation with respect 
to t from the equation 77 = by means of differential consequences of £ and the subsequent 
solving of the resulting equation fj = with respect to its leading derivative u p can be viewed 
as replacing the joint system of £ and r\ = by the equivalent orthonomic system S 

u t = H, u p = f] 

without mixed derivatives on the left hand side. Here the function H coincides with that defined 
after equation ([6]). The leading derivatives of this system are ut and u p ; the principal derivatives 
are u a , where ao ^ 1 or at ^ p; and the other derivatives no, . . . , u p are parametric. 

Proposition 2. The equation £ is conditionally invariant with respect to the operator Q = r\d u 
if and only if the system S is passive with respect to the above ranking. 

Proof. The infinite prolongation of S is the system 5oo 
11 ,1 —n a °n ai H 7; , — D ao D ai f) 

The simplest possibility for deriving integrability conditions of S is to equate the expressions 
for mixed derivatives obtained by differentiating the first and second equations, respectively: 
Dt° Dx 1+P H = -D" " 1 " 1 !?" 1 ^. Each of the derived equations is an identity on equations of 5oo in- 
volving only derivatives lower than the associated mixed derivative (and, consequently, there are 
no other differential consequences) if and only if the conditional invariance criterion is satisfied 
by the equation £ and the operator Q = r]d u , cf. equation ([6]). □ 

5 Relation to involutivity of vector fields 

A connection between generalized conditional symmetries of systems of evolution equations (in 
terms of invariant manifolds) and involutivity of certain system of vector fields was first noted by 
Kaptsov [32] (see also [H p. 131]). For simplicity and uniformity, we restrict our considerations 
to the class (pQ). 
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Let the function u be a solution of the joint system S of the equations £ and u p = fj. We 
introduce the new dependent variables v a ~ l = it a -i and two vector fields 

p-i 

D x = d x + v a d v a-i + fid vP -i , D t = d t + (D a x - 1 H)d v a-i , 

0=1 

where H = H(t, x, v°, . . . , v r ) if p > r, H = H(t, x, v°, . . . , v p ~ l , fj, D x f), . . . , D r x p fj) if p ^ r, 
and n a _i is replaced by v a ~ l in 77. 

In view of the equations for u, the functions u a_1 satisfy the system of differential equations 

v a x - l = v a , a = l,...,p-l, vP- 1 =fi(t,x,v°,...,v' } - 1 ), v b t - l = D b - l H (7) 

which is associated with the system of vector fields {Dt, D x }. 

Proposition 3. The equation E is conditionally invariant with respect to the operator Q = r\d u 
if and only if the system of vector fields {Dt, D x } is in involution. 

Proof. [D t ,D x ] = {D X H — D t fj)d vP -i. Therefore, the system of vector fields {D t ,D x } is in 
involution if and only if these vector fields commute, i.e., D X H — Dtf] = 0. This last equation, 
after the inverse substitution u a _i = v a ~ l , is equivalent to equation ©. □ 

If the system of vector fields {D t , D x } is in involution, the associated system © is completely 
integrable in the old terminology (see, e.g., [TTJ p. 1]). 

Corollary 1. A (1 + 1) -dimensional evolution equation £ is conditionally invariant with respect 
to a pth-order operator Q in reduced form if and only if it possesses a p-parametric family of 
Q -invariant solutions. 

Proof. Suppose that the equation £ is conditionally invariant with respect to the operator Q. 
Then the system of vector fields {Dt,D x } is in involution and, therefore, is integrable by 
the Probenius theorem. The dimension of the span of {Dt,D x } equals two for any fixed 
point (t, x, v°, . . . , v p ~ l ). Therefore, the general solution of the system ([7]) is parameterized 
by2 + /3 — 2 = p arbitrary constants. Its projection to v° necessarily contains all the arbitrary 
constants and gives the general solution of the joint system of £ and fj = 0. 

If the equation £ is not conditionally invariant with respect to the operator Q then the system 
of vector fields {Dt,D x } is not in involution and can be iteratively completed for integrability 
by [Dt,D x ] and the other subsequent commutators which do not lie in the span (over the ring 
of smooth functions) of the system of vector fields from the previous steps. The dimension of 
the span of the completed system is greater than two. (We consider a neighborhood of a point 
in which D X H — Dtf) 7^ 0.) Therefore, the general solution of system (JT]) is parameterized by 
less than p arbitrary constants. □ 

Corollary 2. The set of joint solutions of an equation u p = fj(t,x,U( p _i iX \) and an evolution 
equation £ is parameterized by at most p constants. 



6 Reduction and conditional symmetry 

In this section we discuss ansatzes for the unknown function u, i.e., specific forms for finding 
families of solutions. We shall focus on the following class of (generalized) ansatzes: 

u = F{t,x,Cp{uj)), (p = {ip l ,...,tp p ), (8) 
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where ip 1 , . . . , <p p are new unknown functions of the single invariant variable u = t, det<J> ^ 0. 
By $ and $ we denote the matrices 



d> = = = (jw)» i> = c#° 6 ) = s- 1 . (9) 

Here F a _! = d^F/dx^ 1 and F a _ x ^ b = d a F/dx a ~ 1 d^> . 

Ansatzes u = F l (t, x, <p l {oj)) and u = F (t, x, (p 2 {ui)) with the same number of new unknown 
functions and the same oj = t are called equivalent if there exists a vector-function £ = (^(i,^ 1 ) 
invertible with respect to (p l such that F 2 (t, x, £(t, (p 1 )) = F (t, x,(p l ). This notion of equivalence 
can be extended, e.g., by permitting dependence of (p 1 and (p 2 on different arguments cj% = u)\(t) 
and ll>2 = u>2(t), respectively, but we do not consider this possibility in order to retain the 
distinguished role of the variable t for evolution equations which is fundamental for the general 
line of argument in this paper. 

Lemma 1. Up to the equivalence of ansatzes, for any fixed p there exists a bijection between 
operators of the form ([3|) and ansatzes of the form dH) . 

Proof. An ansatz constructed with an operator Q of the form (J3J) is a representation of the 
general solution of the ordinary differential equation u p = fj(t,x,ur p _i tX )) (with t playing the 
role of a parameter) and, therefore, has the form ([8]) . Equivalent ansatzes only amount to 
different representations of the general solution. (This in fact is the reason for our notion of 
equivalence of ansatzes.) 

The function 77 from the constraint corresponding to an ansatz of the form (JSj) can be 
calculated by the standard method of reconstructing the right hand side of an ordinary dif- 
ferential equation from its general solution. Namely, differentiating the ansatz with respect 
to x up to order p — 1 and solving the resulting system u a _i = F a —i with respect to <p, we 
obtain expressions for (p as a function of t, x and ur p _ ljX y. (p a = T a (t,x,Ur p _i iX \). (This 
is possible since det $ 7^ 0.) Then the ansatz corresponds to the constraint u p = fj, where 
f) = F p (t,x,^l a=Ia{txUip _ hx)) . □ 

The reduction procedure with ansatz ([8]) is implemented in the following way. The substitution 
of dS|) into £ gives the equation F t + F^atpf = H, where H = H(t, x, <p) = H(t, x, Fr TjX \(t, x, (p)). 
We differentiate this equation with respect to x up to order p — 1 and solve the system so 
obtained with respect to (p% (which is possible since det $ 7^ 0). This procedure results in the 
system 

tp* = G a := $ a \H - F t )b-i. (10) 
In general, the right hand sides G a of the equations of this system will be functions of t, x and <p. 

Definition 2. If all the functions G a are independent of x, the system (pf = G a (t, (p) is a well- 
determined system of ordinary differential equations in (p, which is called the reduced system 
associated with the equation £ and ansatz ([8]). In this case we say that the ansatz ([8]) reduces 
the equation £. 

Remark. There also exists another notion of reduction in which a split with respect to the 
independent variables complementary to the invariant ones is possible after substituting ansatzes 
into the initial equations [45J. This kind of reduction is connected with the notion of weak 
symmetry [45, 52] and may be called weak reduction. In contrast to it, Definition [2] gives a special 
case of the general notion of reduction which does not involve a split [431 IUT] . It generalizes the 
classical Lie reduction based on Lie symmetries [121 [46] and the reduction procedures related to 
nonclassical [3 [70] and generalized j27[ section 18.2] symmetries. 
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To allow for a smooth presentation of the subsequent results we now introduce some notions 
related to parametric families of functions and prove some auxiliary statements. 

Definition 3. The parameters >q, . . . , x p are essential in a parametric family {f(t,x, x)} of 
functions of t and x if there do not exist a function / of p + 2 arguments, where p < p, and 
functions £ , . . . , of x such that /(i, x, x) = /(£, x, C 1 (-^ r )> • • • > C p (&))- 

Lemma 2. Zei J 7 = {u = f(t, x, x)} be a parametric family of solutions of £ . All the parameters 
xi, . . . , x p are essential in F (i.e., F is indeed a p-parametric family) if and only if 

d d(/o,...,/ p -l) 
d(xi, . . . , H p ) 



Proof. Suppose to the contrary that all the parameters in T are essential but condition (jlip 
in not satisfied. The latter implies that the values t, x, /o, . . . , f p —\ are functionally de- 
pendent. Thus there exists p' and a function rf of p' + 2 variables such that p' < p and 
/p/ = rj'(t,x, f(p'-i.x))- This means that any solution of £ from the family T also is a solution of 
the equation u p i = rf(t, x, uip)_i x \). Therefore, in view of Corollary [2] the number of essential 
parameters of J- is not greater than p', contradicting our assumption. 

Conversely, if some of the parameters *c\ , . . . , x p are inessential in T then the determinant 
from (jlip must obviously vanish. □ 



Roughly speaking, the parameters in families of solutions of evolution equations are essential 
if and only if they are essential with respect to x. This provides further evidence that (1 + 1)- 
dimensional evolution equations are closely related to ordinary differential equations and in 
various aspects the variable t plays the role of a parameter. 

Definition 4. Families {u = f(t, x, x)} and {u = f(t, x, x')} of functions with the same num- 
ber of parameters are defined to be equivalent if they consist of the same functions and differ 
only by parameterizations, i.e., if there exists an invertible vector-function £ = C(x) such that 
/(i,z,C(x)) = f(t,x,x). 

Now we present the main statements of this section. 

Theorem 1. Up to the re-parametrization equivalence of solution families and the equivalence 
of ansatzes, for any equation of the form (pQ) there exists a one-to-one correspondence between 
p-parametric families of its solutions and ansatzes reducing this equation. 

Proof. Suppose that an ansatz of the form ([8]) reduces £. Since the reduced system is a normal 
system of p first-order ordinary differential equations in ip, its general solution can be represented 
in the form (p = tp(u), x), where x = (xi, . . . , x p ) are arbitrary constants and det(^>* ) 7^ 0. This 
representation is unique up to re-parametrization. Substituting this solution into the ansatz 
results in the p-parametric family T of solutions u = f(t, x, x) of £ with / = F(t,x, ^(t, x)). 
All the parameters >c\ , . . . , x p are essential in T by the chain rule since 

det(/ a _ 1)X6 ) = det(F a „ l5 ^<)| = det$|_^ (ti _ } det%) / 0. (12) 

Conversely, let T = {u = f(t,x, x)} be a /9-parametric family of solutions of £. In view of 
Lemma [2] the expression u = f(t,x,<p(uj)), where uj = t, defines an ansatz for u. This ansatz 
reduces £ to the system (p® = 0. Indeed, after substituting the ansatz into £ we obtain 

(ft + Utf - H(t, f(r, x) ))\- x=m = U\-„ =m <P? = (13) 

since ft = H(t,x, /( r>x ))- We differentiate the last equality in ([TBI) with respect to x up to 
order p — 1 and solve the resulting system with respect to (pt- This system has only the zero 
solution since det(/ _i jX() ) / 0. □ 
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Theorem 2. A (1 + 1) -dimensional evolution equations is conditionally invariant with respect 
to a pth-order evolution vector field Q in reduced form if and only if an ansatz constructed with Q 
reduces the equation £ to a normal system of p first-order ordinary differential equations in the 
p new unknown functions yr, . . . , (p p . 

Proof. Suppose that the equation £ is conditionally invariant with respect to the vector field Q. 
In view of Corollary [IJ the equation £ possesses a p-parametric family {u = f(t,x,x)} of Q- 
invariant solutions. Then the expression u = f(t,x,<p(u)), where u = t, defines an ansatz for u 
associated with Q and reducing the equation £, cf. the proof of Theorem Q3 

Conversely, suppose that an ansatz of the form ([8]) reduces the equation £. Let Q be the oper- 
ator of the form ([3]) associated with this ansatz. Such an operator always exists (cf. Lemma [1]). 
In view of Theorem [1] the ansatz gives a p-parametric family T of joint solutions of the equa- 
tions £ and Q. Then Corollary Q] implies that the equation £ is conditionally invariant with 
respect to the operator Q. □ 

Corollary 3. Up to the re-parametrization equivalence of solution families, for any equation 
of the form ([I]) there exists a one-to-one correspondence between p-parametric families of its 
solutions and canonical pth-order conditional symmetry operators. Namely, each operator of 
this kind corresponds to the family of solutions which are invariant with respect to this operator. 
The problems of the construction of all p-parametric solution families of equation ([I]) and the 
exhaustive description of its canonical pth-order conditional symmetry operators are completely 
equivalent. 

Proof. It is enough to combine Theorems [1] and [2j Each solution constructed with an ansatz 
of the form ([8]) is invariant with respect to the canonical pth-order reduction operator of £ 
associated with the ansatz. □ 

Example 1. Analyzing the results from [28J on the group classification of (1 + l)-dimensional 
variable-coefficient nonlinear diffusion-convection equations of the general form 

f{x)u t = {g(x)A{u)u x ) x + h(x)B(u)u x , 

where f(x)g(x)A(u) ^ 0, we obtain only one essentially variable-coefficient equation 

x 2 u t = (u~ 6/5 u x ) x + x 2 u x (14) 

which is invariant with respect to a realization of the algebra sl(2,R). All the other sl(2,R)- 
invariant equations from the class under consideration are similar (i.e., mapped by point trans- 
formations) to the well-known ("constant-coefficient") Burgers and u _4//3 -diffusion equations. 
Instead of equation (fl~4l) it is more convenient to study the equation 



5 

x 2 v t = vv xx - -(v x ) 2 + x 2 v x (15) 



for the function v = ti -6//5 , i.e., u = i> _5//6 . The maximal Lie invariance algebra of equation (|15p 
is g = (dt, tdt+xd x + 3vd v , t 2 dt + (2tx+x 2 )d x + 6(t + x)vd v ). Extending Lie ansatzes constructed 
by one-dimensional subalgebras of q, we derive the generalized ansatz 

v = 2x 3 + ^{t)x A + ip 5 (t)x 5 + /(t)x 6 , (16) 

which reduces equation (fT5|) to the system of ordinary differential equations 



4 = 7<p 5 - ^ 4 ) 2 , ^ = 18/ - j^V, ^ = - jV) 2 + 2*V- 
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This ansatz represents the general solutions of the equation 

x 3 v xxx - I2x 2 v xx + 60xv x - 120v + I2x 3 = 0. 

In view of Theorem [2j the reduction of equation (115p with the ansatz (|16p is equivalent to the 
fact that (|15p is conditionally invariant with respect to the third-order evolution vector field 

(x 3 v xxx - I2x 2 v xx + 60xv x - 120v + 12x 3 )d v . 

7 No-go theorem on determining equation 

In terms of local solutions, Corollary [3] means that there exists a (local) one-to-one correspon- 
dence between solutions of the determining equation ([6]) and p-parametric families of solutions 
of the initial equation (pQ). We show that this correspondence is realized by transformations 
between systems associated with these equations. 

Theorem 3. The system in the functions 9 a = 6 a (t, x, uq, . . . , u p —i), which consists of the partial 
differential equation ([6]), where fj is identified with 9 P , and the algebraic equations 6 1 = u%, . . . , 
6 p ~ l = Up— i, is reduced by the composition of the differential substitution 

e = (17) 

where X = I(t, x,Uq, . . . , u p -\) are the new unknown functions, := (I® ) and det ^ ^ 0, and 
the hodograph transformation 

the new independent variables: t = t, x = x, x a = I a , 
the new dependent variable: v^ 1 = 

to the system formed by the initial equation £ in the function u = u(t,x,>c) and the equations 
v^ 1 = d b ~ 1 u/dx b ~ 1 , b = 2, . . . , p, where >c plays the role of a parameter tuple and u is identified 
with v° . 

Proof. At first we construct a direct transformation. Extending equation ([6]), we introduce the 
notation 6 1 = u\, . . . , 9 P ~ 1 = u„-%, 9 P = f]. This notation is natural since in view of the 
definition of the operators Dt and D x and equation (|6|), the functions 9 a satisfy the conditions 
DtO a = D X H. We consider the system consisting of the partial differential equation (|6|) and the 
algebraic equations 6 1 = m, . . . , 6 p ~ l = u p _i and carry out the differential substitution (fT7|) . 
In other words, I is a tuple of solutions of the equation D X X = with det(X" ) / 0. It is 
determined by up to the transformation X — > G(t,I), where Gj b ^ 0. Then we carry out the 
hodograph transformation (|18j) . In what follows, for convenience we denote the function v° by u 
and the derivatives d k u/dx k by Uk, k = 1, 2, ... . Differentiating the equality x = X with respect 
to x, we obtain Z x + v b ~ 1 I Ub _ 1 = 0. As D X Z = and det^/ ^ 0, this means that v- -1 = v b , 
b < p, v ( ~ 1 = f](t, x, v°, . . . , v p ~ v ), and therefore v b " 1 = Ub-i an d u p = fj = rj(t, x, Ur p -i tX \), i.e., 
9 a = u a . In the new variables we also have that D x = d x + (D x I a )d Ha = d x . This operator acts 
on the functions of t, x and derivatives of u as the operator D x of total derivation with respect 
to the variable x. Hence D x fj = u p+ k and H = H := H(t,x,ur r)X )). Analogously 

A = d~ t + {D t l a )d Ha = d- t - {u b _ lf - Dt x E)X a Uh _^ a 

since Xf + tt fc _ 1 fZ^. = 0. Moreover, as u p = fj, we also have u pKa = D^ a fj = fju b _ 1 Ub_i^ a , and 
the matrix (ub—i,x a ) ^ s ^ ne i nverse of the matrix (1^ ). This is why in the new variables the 
equation ([6]) takes the form 

De(u i -H) = fjz b _ 1 Dl- 1 (u i -H). 
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For a fixed function u, the equation w p = rju b _ 1 Wb-i with respect to the function w = w(t, x, x) is 
a pth-order ordinary differential equation, with x as the independent variable and t and x playing 
the role of parameters. The functions u Ka are linearly independent solutions of this equation 
since det(£t Xa) {,_i) 7^ 0. Therefore, there exist functions Q a = ( a (t, x) such that tt f - — H = C^u^. 
In view of the indeterminacy of X, we can make the transformation x — >■ x) to transform 
the last equation to the equation of the same form with £ a = 0. 

Conversely, let u = u(t, x, x) be a p-parametric solution of equation ([1]). (We use the notation 
with tildes to be consistent with the first part of the proof.) Assuming v b ~ l = Ub~i as the un- 
known functions, we obtain the system t;~ = H(t, x, v9 -■,) and v ~ _1 = v b , b < p. We successively 
carry out the inverse of the hodograph transformation (|18p and the inverse of the differential 
substitution (fT71) and denote the function 9 P = 9 p (t,x,uo, . . . , u p _i) by fj. By construction we 
have that 6 l = u\, . . . , O^ 1 = u p -\ and for each x the solution u = u(t, x, x) of ([1]) is invariant 
with respect to the operator Q = {u p — fj)d u . This means that Q is an operator of generalized 
conditional symmetry of ([T]) and, therefore, the function fj satisfies equation ([6]). □ 

We call Theorem [3] "a no-go theorem" since it basically states that solving the determining 
equation for generalized conditional symmetry operators is as difficult as solving the original 
equation. It generalizes the analogous no-go theorem on the determining equations for usual 
conditional symmetry operators of evolution equations, whose coefficient of dt is equal to zero 
[17\ [Ml HHl H9| l62l 69j . The main problem in generalizing that result was that the corresponding 
hodograph transformation should involve p independent variables. At the same time, both the 
initial and determining equations involve only a single dependent variable. 

Note that the attribute "no-go" should be treated as impossibility of the exhaustive solution of 
the problem. At the same time, imposing additional constraints on the differential function i) = 
f](t, x, ii(p_i : x)) or choosing a specific form for this function, one can construct a number of 
particular examples of generalized conditional symmetries and then apply them to finding exact 
solutions of the original equation £. Since the determining equation ([6]) has more independent 
variables and, therefore, more degrees of freedom, often it is more convenient to guess a simple 
solution or a simple ansatz for the determining equation, which may then provide a parametric 
set of more complicated solutions of the original equation £. 

This situation is similar to that of Lie symmetries of first-order ordinary differential equations. 
Indeed, the solution of the determining equation for Lie symmetries of a first-order ordinary 
differential equation C is a much more complicated problem than the solution of the original 
equation C Even if a Lie symmetry generator of C is known, it may be just as difficult to find 
an invariant of the associated one-parameter group (which is a necessary step of solving by the 
Lie method) as it was to integrate the original differential equation C [421 pp. 131-133]. At the 
same time, certain first-order ordinary differential equations (e.g., homogeneous ones) possess 
simple Lie symmetries which can easily be found by an educated guess and then effectively used 
for the integration of these equations. 

The above approach to the construction of exact solutions using generalized conditional 
symmetries of special kinds was applied in the literature to a number of different classes of 
evolution equations, in particular to quasilinear second-order evolution equations. We recall 
only some of these results. 

Generalized conditional symmetries of many particular cases of equations of the general form 

ut = g 3 (t,x,u)u xx + g 2 (t,x,u)u 2 x + g x (t,x,u)u x +g°(t,x,u) 

were looked for by a number of authors in a form similar to the right hand sides of the corre- 
sponding equations, 

r] = u xx + g 2 (t,x,u)u 2 x +g 1 (t,x,u)u x +g°(t,x,u), 
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or in the equivalent form rj = n t +5 2 (t,x,ti)n^ + 5 1 (t,x,u)ti x + 5°(t,2;,u), see, e.g., [22 | \30 \ [29 } |53" 1 
[54"1 [67] and references therein. Another intensively investigated class of generalized conditional 
symmetries and generalized ansatzes is related to differential constraints which are equivalent 
to linear differential constraints with respect to point transformations, see, e.g., [9j 120 1 [25 | [67] 
and references therein and cf. also Example [TJ 

As shown in the next section, a generalized first-order conditional symmetry in canonical 
form (u x — rj(t, x, u))d u of an evolution equation is, up to sign, the evolution form of the singular 
nonclassical symmetry operator d x + r/(t,x,u)d u of the same equation. In |51] such symmetries 
of different classes of quasilinear second-order evolution equations were studied under the as- 
sumption of separation of variables in the coefficient rj, n = ^°(t)^ 1 (x)C 2 (w). Earlier the partial 
case rj = C 1 (x)C 2 («) was investigated in [23] for equations of the form it t = u%u xx + p>u J +1 + f(u). 
The important special subcases C 1 ( x ) = x an d C l { x ) = x_1 were separated therein. The lat- 
ter subcase, which generalizes scale-invariant solutions, was considered within a more general 
framework in |21j . An extension of results obtained in |23] was presented in [55]. The ansatz 
f] = r] l {t,x)u a+1 + 7]°(t,x)u a was used in [26] for the fast diffusion equations of the form 
u t = (u~ a u x ) x . 

8 Usual and generalized reduction operators 

It seems natural that usual conditional symmetry is a particular case of generalized conditional 
symmetry. On the other hand, the criterion of usual conditional symmetry restricted to the 
case of evolution equations is essentially different from (j4| . This is why we formulate the precise 
relation as a proposition. 

Proposition 4. The vector field Q = rdt + t,d x + r]d u , where the coefficients t, ^ and r] are 
functions oft, x and u, is a usual conditional symmetry operator of an equation £ of the form §\§ 
if and only if the operator Q = i)d u , where r) = t\—tB.—^u x , is a generalized conditional symmetry 
operator of the same equation. 

Proof. The first way of proving this is simpler but essentially involves statements on properties 
of the corresponding families of invariant solutions. A solution of £ is Q-invariant if and only if 
it is Q-invariant. Moreover, ordr) = p, where p = r if r 7^ and p = 1 if r = 0. Suppose that 
Q is a usual reduction operator of £. Propositions 2 and 5 from [3l] imply that the equation £ 
possesses an r-parametric (resp. one-parametric) family of Q-invariant solutions if r 7^ (resp. 
r = 0). Then Corollary Q] implies that Q is a generalized conditional symmetry operator of £. 
The proof of the converse is similar. 

The second way is more direct and technical. We have to show that the corresponding 
invariance criteria are equivalent. In what follows E = ut — H ', rj = Q[u] = 77 — rut — £,u x , 
k = 0, . . . , r, k = 0, . . . , p, j = 0, . . . , r — 1 and j = 0, . . . , p— 1. We have {Qt r \ — Qt r ))E = ED t r — 
£,D X E — H Uk D^E . The expression Q( T )E involves at most the derivatives Uk and u t Hence the 
differential consequences which should be taken into account in the usual conditional invariance 
criterion are exhausted by £ itself and D 3 x fj = 0. Analogously, the expression Q^(u t — H) 
involves at most the derivatives u m , m = 0, . . . , r + p, and u t ^. Therefore, the differential 
consequences which should be taken into account in the usual conditional invariance criterion 
are exhausted by D^E = and D^r) = 0. Finally, we have the chain of equivalences 




£r + p{~}Q.(r+p,x) 



= 



and the result follows. 



□ 
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Despite the fact that the sets of Q- and Q-invariant solutions of £ coincide, in the case t / 
the procedures of the reduction of £ with respect to the operators Q and Q are quite different. 
Thus, the operator Q reduces £ to a single rth-order ODE in a single unknown function, where 
the invariant independent variable necessarily depends on x or u. The operator Q reduces £ to 
a system of r first-order ODE in r unknown functions, where t can be taken as the invariant 
independent variable. We illustrate this situation by the following example. 

Example 2. Usual reduction operators of variable coefficient semilinear diffusion equations with 
power sources were investigated in [63, 64]. Namely, the equations studied have the general form 



f{x)u t = {g{x)u x ) x + h{x)u m , (19) 

where /, g and h are arbitrary smooth functions of the variable x, f(x)g(x)h(x) ^ 0, and m is 
an arbitrary constant, m/0,1. The most convenient approach to this problem, as it turns out, 
is to map the class (fl~9j) to the class 



v t = v xx + H{x)v m + F(x)v (20) 

by a family of point transformations parameterized by arbitrary elements /, g and h and then to 
investigate usual reduction operators of equations from the latter class. Under both the group 
classification and the classification of reduction operators the following equation is singled out 
from the class (1201): 



v 3 9 v 

Vt = Vxx ~x^ + l^' (21) 



Note that by the point transformation t = t, x = In \x\, v = vj \j2\x\ equation (f2Tj) is reduced 
to the equation e 2x vj: = v xx — 2v 3 + 2v. 

The maximal Lie invariance algebra of (|2ip is generated by the operators D = Atdt+2xd x +vd v 
and dt- Inequivalent non-Lie usual reduction operators of (|2ip . having nonzero coefficients of dt 
are exhausted, up to the discrete symmetry transformation of alternating the sign of v, by the 
operator 



Q = d t + 



v 3 \ „ 3 / v 3\/2 v 2 v 

~x) x ~2\x^~ ~l~x^ ~x^ + ^J 2 ) v ' 



For all expressions to be well-defined, we have to restrict ourselves to values x > 0. (Another 
way is to replace x by \x\.) 

We discuss two ways of using the operator Q for finding exact solutions of equation ([2T|) . 

First way. To construct an ansatz with the operator Q, we have to solve the quasilinear first- 
order partial differential equation Q[v] = 0. The corresponding invariant independent variable 
necessarily involves the dependent variable v. For simplifying calculations, we suppose at first 
that v t ^ and carry out the hodograph transformation t = v, x = x, v = t which maps 
equation ([2"T]) and the reduction operator Q to the equation 

t 3 9 [ 

Vt 2 v xx + v x 2 v u - 2 v~ t v x v~ t - + v- 2 + -Tg - - -rj Vt 3 = ( 22 ) 
and its reduction operator 

A _ 3fP 3y/2 P i , 2y/2\ a /^3\/2^_3\ 

V- 2VS 3 2 x5/2 5-2 + ^p+y 2 ^3/2 zJ^ + Ov, 
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respectively. An ansatz constructed with the operator Q has the form 

, . 1 „ot + V2x~ 1 _ 2 , „ 2 i-y/2x 
v = z{oj) H x == x , where ui = x ==, 

and reduces ([22]) to the single second-order linear ordinary differential equation a;,8 ww + 2z w = 
in the function z = z(ui). After substituting to the ansatz, the general solution z = c\ + C2UJ~ 1 
of the reduced equation gives the exact solution 

_ x 4 + 24c 2 t + y/Zt 1 n ~ 



24x 2 t-VM 12* 

of (|22p . Applying the inverse hodograph transformation, we construct the non-Lie solution 

, = v^ 3 - 4 4 + / 24t + C f 2 2 - C2 (23) 
x 4 + (24t + ci)x 2 + c 2 V ; 

of equation (|21|) . where c\ = — 24ci and C2 = — 24c2- The constant c\ can be canceled by a 
translation with respect to t. If C2 7^ 0, this constant can be set to 1 by a scale transformation 
generated by D. (Recall that the above transformations are Lie symmetries of equation (|2ip .) 
The solution (|23p with C2 = is a Lie solution invariant with respect to the scale symmetry 
operator D. However, it is much harder to find this solution by the reduction with respect to 
the operator D. The corresponding ansatz v = \^2xz(oj), where u = x/y/\t\, has a simple form 
but the reduced ODE ^z^ + lo(1 — uS)z w + 2z — 2z 3 = is nonlinear and complicated. 

Under the condition vt = the equation Q[u] = implies equation (I21D and is reduced by 
the transformation v = \/2xz(x) to the equation (z — l)(xz x + z 2 — 1) = which is equivalent 
to the Riccati equation xz x = 1 — z 2 . Therefore, the corresponding invariant solutions of (|21|) 
have the form 

v = ^Tx ClX \~ C \ (24) 

where only the ratio of the constants c\ and C2 is essential. Note that a function v is a stationary 
solution of (|2ip if and only if v = y2xz(x), where z = z(x) is a solution of the equation 
z xx = 2(z 3 — z) which is integrable in terms of elliptic functions. 

Second way. Another way to use the operator Q for finding exact solutions of equation (|21l) 
is to consider the second-order reduction operator rj(t,x,v,v x ,v xx )d v , where the differential 
function coincides with the characteristic Q[v] on the manifold determined by equation (I2ip in 
the corresponding second-order jet space. Here 

3^2 vv x 3 1 v 3 9V2 v 2 3 v 3^/2 
V - -Vxx - —^/2 + ~ Vx ~ 2^3 + ~^~^5/2~ ~ 4^2 ~ ^372' 

The associated invariant surface condition rj = is a second-order ordinary differential equation, 
where x and v are independent and dependent variables, respectively, and t plays the role of a 
parameter. It is reduced by the differential substitution 



to the linear equation 

x 3 w xxx - 3xw x + 3w = (25) 
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in the new unknown function w = w(t,x), whose general solution is w = ip^^x 3 + ip 1 (t)x + 
ip 2 (t)x . Therefore, we have the following ansatz for the function v: 

v _.^ W°(t)x* + 4>Ht)x 2 -nt) (26) 

where only two ratios of the functions ip^, fj, = 1, 2, 3, are essential. 

To make a conventional reduction of equation (|2ip with ansatz (|26p . we would suppose that 
one of the functions tp 11 , e.g. ifi°, is nonvanishing. After substituting ansatz (|26|) into (|2ip . we 
would obtain the reduced system of two first-order ODEs in the functions (p l = ip l /ip® , i = 1, 2. 
Then it would be necessary to consider the case ip° = and ip 1 ^ and to derive the reduced 
first-order ODE in tp = ip 2 /ip l . The condition (ip , ip 1 ) = (0,0) leads to the single solution 
v = —y/2x. This partition into different cases corresponds to the partition made in the first way. 

We use a more advanced technique allowing us to avoid the consideration of different cases. 
The entire systems of the equation r\ = and equation (|2Tj) is equivalent to the system of the 
equations ([25]) and 



w t = 3w xx + \w x - ^w. (27) 



x 



Moreover, (x 3 w xxx — 3xw x + 3w)d w is an operator of generalized conditional symmetry of (f2"T|) . 
Therefore, the associated ansatz w = ip°(t)x 3 + tp 1 (t)x + ip 2 (t)x~ 1 reduces equation (|27l) . and the 
corresponding reduced system is ip® = 0, ip\ = 24^°, ip 2 = with the general solution ip° = Co, 
tp 1 = 24cQt + ci, ip 2 = C2- As a result, we have the solution 



v 



3c x 4 + (24c t + d)x 2 - c 2 

X cqx 4 + (24c i + ci)x 2 + c 2 



of equation ([2?]) . The conditions Co 7^ and cq = correspond to the solutions ([23]) and ([24"|) . 
respectively. 

In a similar way, the conversion of usual nonclassical symmetries into generalized ones was 
implicitly used, e.g., in [21 [101 E] in the reduction of the nonlinear (constant coefficient) heat 
equations with source terms in the form of cubic polynomials, including the Fitzhugh-Nagumo 
equation. 



9 Conclusion 

In the study of generalized conditional symmetry of evolution equations we adhere to the fol- 
lowing principles: 

1. The property of an operator Q to be a conditional symmetry of a differential equation C is 
equivalent to the fact that the corresponding invariant surface equation Q[u] =0 is formally 
compatible (in a certain sense) with C, i.e., the joint system of the above two equations has no 
nontrivial differential consequences. This determines what differential consequences of these 
equations should be involved in the criterion of the conditional invariance of the equation C 
with respect to the operator Q. In the property of formal compatibility conditional symme- 
tries differ from purely weak symmetries [45, 52J for which the calculation of integrability 
conditions (resp. "the reduction to passive form") of the corresponding joint systems has to 
be carried out in each case. In fact, weak symmetries of C are associated with differential 
constraints whose solution sets at least intersect the solution set of C 
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2. Therefore, the criterion of conditional invariance in fact is nothing but the criterion of formal 
compatibility for a system associated with the pair of equations C and Q[u] = 0. This has 
two consequences: There does not exist a universal explicit criterion of conditional invariance 
similar to the criterion of Lie invariance, which would contain a priori the complete informa- 
tion which differential consequences to take into account and would be appropriate for any 
system of differential equations and any set of generalized vector fields. At the same time, 
for any fixed pair of a system of differential equations and a set of generalized vector fields 
the criterion can be formulated in different forms. 

3. Single generalized conditional symmetries are assumed equivalent if they differ by multipliers 
being nonvanishing differential functions. Therefore it suffices to consider only symmetries 
with characteristic containing some isolated (e.g., highest-order) derivative of the unknown 
function. 

4. In order to be usable, a conditional symmetry should correspond to an integrable differential 
constraint which admits a simple representation of its general solution. Such a representation 
is considered as an ansatz for the solution of the initial equation C. The formal compatibility 
of the differential constraint with C should imply a (strong) reduction of C by the ansatz. 
In other words, after the substitution of the ansatz into C we should obtain a system of 
differential equations of a simpler structure, e.g., with a smaller number of independent vari- 
ables. Symmetries equivalent as vector fields induce the same set of ansatzes and equivalent 
reductions. In fact there does not exist a universal precise definition of reduction which does 
not involve splitting with respect to parametric variables and covers all possible represen- 
tations of solutions. In view of the above problems of integrability and defining reduction, 
it is still unclear in the general case what differential constraints formally compatible with 
the initial equation should be considered associated with reduction operators. This question 
becomes trivial and has a positive answer in the situation considered in the paper (single 
evolution equations and differential constraints depending only on derivatives with respect 
to x). Probably, in the general case it would be more natural to assume that the notion of 
reduction operator is narrower than the notion of conditional symmetry, cf. [43j. 

5. If the characteristics of operators coincide on the manifold determined by the initial equa- 
tion C and one of the operators corresponds to a differential constraint formally compatible 
with C then the other operators have the same property. Such conditional symmetry oper- 
ators can be considered equivalent in a weak sense since they are associated with the same 
set of invariant solutions of C. At the same time, they are inequivalent, in general, from the 
point of view of their usefulness for finding solutions. In particular, they may give inequivalent 
ansatzes and reduced systems. 

These general principles can be applied in other situations as well. We plan to complete soon 
our study on basic properties of usual (i.e., first-order quasilinear) conditional symmetries of 
systems of differential equations. 

In spite of the no-go results presented in the paper, generalized conditional symmetries can be 
effectively applied to the construction of exact solutions of evolution equations. As it is impos- 
sible to exhaustively describe generalized conditional symmetries of a fixed evolution equation, 
they should be looked for under additional constraints or in special classes of differential func- 
tions, e.g., with separated variables. In this way, usual and generalized conditional symmetries 
were studied for a number of particular subclasses of evolution equations, cf. the discussion in 
the end of Section [7J Note that only in [33] generalized conditional symmetries which are not in 
reduced form were considered. Generalized conditional symmetries were also used for the exact 
solution of initial- value problems for evolution equations [U [68] . Another relevant direction of 
research is the related inverse problem, namely, the description of evolution equations possessing 
certain generalized conditional symmetries, see [311 l58l l59j [BT] and references therein. 
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A systematic investigation of generalized conditional symmetries of non-evolution equations 
in fact is not available in the literature at the moment. An exception is the paper |43j of Olver 
mentioned in the introduction, where the connection between the reduction of a partial differen- 
tial equation by a generalized ansatz within the higher-order direct method of Galaktionov [20] 
and the compatibility of the associated differential constraint with this equation was discovered. 
At the same time, there exist a number of examples on the application of generalized ansatzes 
to finding exact solutions of non-evolution equations, which are collected, e.g., in [151 125] . It 
is obvious that all such examples can be interpreted within the framework proposed in [43] . 
Ansatzes of another kind with new unknown functions depending on different arguments arise 
under generalized separation of variables [U [66]. Theoretical aspects of this subject should 
certainly be further investigated. 
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